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^ . ABSTRACT 

' Dust plays various important roles in galaxy formation and evolution. In the early epoch of 

galaxy evolution, dust is only supplied by supe rnovae (SNe). With the aid of a new physical 
7— I ■ model of dust production by SNe developed by Nozawa et al. (2003)(N03), we constructed a 

J> ' model of dust emis sion from forming galaxies on the basis of the theoretical framework of 

, iTakeuchi et alJ ( 120031) (T03). N03 showed that the produced dust species depends strongly on 

the mixing within SNe. We treated both unmixed and mixed cases and calculated the infrared 
(IR) spectral energy distribution (SED) of forming galaxies for both cases. Our model SED 
is less luminous than the SED of T03 model by a factor of 2-3. The difference is due to our 
improved treatment of UV photon absorption cross section, as well as different grain size 
| and species newly adopted in this work. The SED for the unmixed case is found to have 

-»■»« . an enhanced near to mid-IR (N-MIR) continuum radiation in its early phase of the evolution 

' (age < 10 7 25 yr) compared with that for the mixed case. The strong N-MIR continuum is due 

to the emission from Si grains, which only exist in the species of the unmixed dust production. 
Q ■ We also calculated the IR extinction curves for forming galaxies. N03 dust was found to yield 

| a smaller extinction than that of T03 model. For the unmixed case, near-IR (NIR) extinction is 

r/3 , dominated by large grains of Si and amorphous carbon, and silicate features are less prominent 

compared to the curve given by T03. On the contrary, the extinction curve of the mixed case 
£>. | has a similar shape with that of T03. Then we calculated the SED of a local starbursting 

• *h . dwarf galaxy SBS 0335—052. Our present model SED naturally reproduced the strong N- 

y\ ' MIR continuum and the lack of cold FIR emission of SBS 0335—052. We found that only 

\ the SED of unmixed case can reproduce the NIR continuum of this galaxy. We then made a 

prediction for the SED of another typical star-forming dwarf, I Zw 18. The MIR continuum of 
I Zw 18 is expected to be much weaker than that of T03 SED. We also presented the evolution 
of the SED of LBGs. Finally, we discussed the possibility of observing forming galaxies at 
: ..J 

Key words: dust, extinction - galaxies: dwarf - galaxies: ISM - galaxy formation - infrared: 
galaxies 
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1 INTRODUCTION 

In spite of the recent vast progress both in observational and theo- 
retical studies, our understanding of the physics of galaxy forma- 
tion and evolution is still far from sufficient. 

The cosmic star formation (SF) history, introduced by 
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ITinslev & Danlv (1980) and developed by subsequent studies (e.g., 
.Lilly et alJl996tlMadau et alll996t) . has always dra wn much atten- 
tion. Now the observations re ach up to ~ 6-8 (e.g. lStanwav et all 
2003; Bou wens et all2004alo) . In such studies, the role of dust has 
been increasingly recognized when we try to understand the evo- 
lution of galaxies in the context of cosmic star formation history, 
because dust grains absorb stellar light and re-emit it in the far in- 
frared (FIR). Even a small amount of dust can lead to a signifi- 
cant underestimation of the star fo rmation rate (SFR) dSteidel et alJ 
1999; Adelbereer & Steidel 2000). Indeed, there is another extreme 
category of high-z galaxies which have large amount of dust and 
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are extremely lum inous in the FIR an d submillimetre (submm) 
wavelengths (e.g.. hughes et alj[l99ll lEales et al, 120031). Heav- 
ily hidden SF is suggested in the se galaxies (e.g., Takeuc hi et alJ 
l2001albtlTo"tani & Takeuchil2002t) . 

Further, from physical point of view, dust grains are one of the 
fundamental ingredients in the activity of galaxies. Since they are 
formed in a variety of environments ranging from explosive ejecta 
of novae and supernovae to the outflowing gas of evolved low-mass 
stars, the dust formation is closely related to the SF activity iDwekl 
1998). Moreover, the existence of dust is crucial in the physical 
process of galaxy formation and evolution through the formation 
of molecular hydrogen (e.g.. iHirashita & Ferrarall2002h . The dust 
itself plays a leading part in the physics of SF activity in galaxies 
in the whole period of the cosmic history. 

Then, how about very young galaxies in the early universe 
at z > 5? It is often assumed, without deliberation, that the ef- 
fect of dust is negligible for such young galaxies, because of their 
low metallicities. There is, however, a good counter-intuitive ex- 
ample in the local Universe: a local dwarf star-forming galaxy 
SBS 0335—052 has a very young stellar age (< 10 7 yr) and low 
metallicity (1/41 Z@), but has a heavily embedded active star 
formation and strong continuum radiation in N-MIR wavelength 
regime iDaleetal. 2001; Hunt et al. 2001). In addition, there are 
increasing number of observations suggesting the exi stence of dust 
in high-z s ystems such as Lyman a systems (e.g., iLedoux et alJ 
120021120031) or QSOs (e.g.jBertoldi e t all2003l) . 

To produce dust effectively in such a young system, the dust 
enrichment should have occurred primarily in the ejecta of super- 
novae (SNe), especially Type II supernova (SN II) explosions , be- 
cause t he lifetim e of the progenitor is short enough (~ 10 6 yr) (e.g., 
iDwek & Scalolll98CtlKozasa et allll989h . Recent theoretical stud- 
ies claim the possibility of very massive stars for the first generation 
stars, which end their lives as pa ir-instability supernovae (PISNe) 
I Hege r & Woosle\l2002tlUmeda & Nomotol2002h . We should also 
take into account for the PISNe to study the very early evolution of 
dust. Recently, the dust formation in SNe ejec ta has been observa- 
tionally supported (e.g.. lDouvion et aTll200ll : lMorgan et alJl2003l 
iDunne et all2003l) . 

For investigating the properties of dust in details, theoretical 
predictions for the amount and composition of dust are required. 
So far, some theoretical models o f dust production by SNe have 
been developed. Todini & Ferrara 1 200 1 , hereafter TF01) showed 
that the dust mass produced by a SN II is 0.1-0.4 Mp, ap plying 
the the ory of nucleation and grain growth by Ko zasa & Hasegawa 
ll987l) . In the calculations, they assumed an adiabatic cooling in 
the ejecta with uniform gas density as well as elemental compo- 
sition within He-core. They also found that SNe form amorphous 
ca rbon with size around 300 A an d silicate grains around 10-20 
A. ISchneider. Ferrara & Salvaterral i200 4) extended the progenitor 
mass range to the regime of PISNe (140-260 Mq) and found that 
10-60 Mq of dust forms per PISN. The grain radius depends on 
the species and is distributed from 0.001 to 0.3 /J,m. 

Recently, iNozawa et alJ J2003T) (hereafter N03) carefully took 
into account the radial density profile and the temperature evolution 
in the calculation of the dust mass in the ejecta of SNe II and PISNe. 
N03 showed that the produced dust species depends strongly on the 
mixing within SNe. In particular, carbon dust is not produced in the 
mixed case, because the carbon and oxygen are mixed and com- 
bined to form CO molecules. On the contrary, it forms in unmixed 
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SN, since there is a carbon-rich region at a certain location in the 
ejecta of SNe. In addition, N03 predicts a dust mass larger than that 
of TF01 for SNe II. 

The amoun t of dust produced in SNe is thus still a matter 
of debate. Morgan & Edmunds (2003) adopted a dust production 
smaller than those of TF01 or N03 in their theoretical calcula- 
tion, in line with IR observations at the time. On the other hand, 
various observational efforts have been devoted to constrain the 
dust p roduction effi ciency of SNe (e.g..lArendt. Dwek. & Moseley l 
I l999t IDunne et aljEool iMorgan et alJl2003t iGreen et al J 12004 
Mines et all |2004 IKrause et all |2004 IWilson & Batrlal hooH) . 
These results range from 0.003 Mq to 1.0 M© per SN, and N03 
dust production seems consistent with the upper range of these 
observational constraints. However, in these observational evalu- 
ations, a large uncertainty still remains due to the uncertainties in 
the estimation of dust mass in Cas A caused by the foreground con- 
tamination iKrause et al]l2004 . although the evidence for dust in 
some SN remnants (SN1987A and Kepler) is still pointing towards 
SN origin. Thus, further close interaction between theoretical and 
observational works is required to overcome the uncertainties and 
obtain a reasonable picture of dust production in SNe. 

In order to examine the dust properties of such high-z galax- 
ies, the most direct observable is their spectral energy distribution 
(SED), especially at the FIR wavelengths. IHirashita et aD 1 2002 ) 
modeled the evolution of FIR luminosity and dust temperature in 
a young starburs t on the basis of SNe II grain formation model 
of TF01 . ITakeuchi et alJ j2003) (T03) subsequently constructed a 
model of infrared (IR) spectral energy distribution (SED) of galax- 
ies starting from the model of Hirashi ta et alJ I2002I) . T03, for the 
first time, properly consider the dust size distribution peculiar to 
the very early stage of galaxy evolution in the model of the IR 
SED of very young galaxies, and successfully reproduced the pe- 
culiar MIR SED of SBS 0335-052 (tho ugh their result seems to 
have overestimated the FIR continuum). ITakeuchi & Ishiil 12004 
(T04) applied the T03 model to the Lyman-break galaxies. They 
found that the suggested hot dust in the Lyman-break galaxies (e.g., 
lOuchi et alll999MChapman et all20o51lSawickil200ll) can be nat- 
urally explained by the T03 model. 

Since these works are based on TF01 conjecture, the dust 
formation model used is based on the classical nucleation theory 
(Feder et al. 1966). Recently, however, the importance of nonequi- 
librium (non-steady state) effects on the dust grain for mation has 
been recognized in variou s astrophysical contexts (e.g., iGail et all 
ll984lTanaka et all2002t N03). In addition, N03 found that the ra- 
dial density profile of the SN progenitor and the temperature evolu- 
tion of the ejecta also affect the dust grain formation. Hence, now 
is the time to take into account these concepts to further investi- 
gate the dust emission from young galaxies. In this work, we con- 
struct a new model of IR SED of extremely young galaxies based 
on N03 SN dust formation model. Starting from the size distribu- 
tion and the amount of dust predicted by N03, we calculate the dust 
emission model by extending the T03 model to treat multiple dust 
species. 

The paper is organized as follows: In ij2|we explain the frame- 
work of our SED model. We present the basic result, the evolution 
of the SED of extremely young galaxies, in J3| Related discussions 
on local star-forming dwarf galaxies and high-z galaxies will be 
in 21 !0i s devoted to our conclusions. Throughout this paper, we 
use a cosmological parameter set of (h, Qo, Ao) = (0.7, 0.3, 0.7), 
where h = .Ho/100 [kms _1 Mpc -1 ]. 
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Figure 1. Dust size distributions of type II supernova model with the progenitor mass of 20 Mq bv lNozawa et alJ 120031) . The thick lines are the total dust 
size distribution function. Labels 'Unmixed' and 'Mixed' indicate the adopted assumption of mixing of elements in the ejecta. We assumed 20 Mq as a 
representative for the supernova progenitor mass in this work. 



Table 1. Dust grain species. 



Species 


Unmixed 


Mixed 


Density (p;) 
[gem -3 ] 


Ref" 


C 







2.28 


1 


Si 







2.34 


2 


Si0 2 








2.66 


3 


Fe 







7.95 


4 


FeS 







4.87 


5 


Fe 3 4 







5.25 


6 


AI2O3 








4.01 


7 


MgO 







3.59 


8 


MgSi0 3 








3.20 


9 


Mg 2 Si0 4 








3.23 


5 
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2 SED MODEL FOR FORMING GALAXIES 

2.1 Species and size distribution of dust grains produced by 
SNell 

2.1.1 Dust production model ofNozawa et al. (2003) (N03) 

N03 investigated the formation of dust grains in the ejecta of Pop- 
ulation III SNe (SNe II and PISNe, whose progenitors are initially 
metal-free). As we mentioned above, they treat some aspects which 
TF01 have not taken into account: (i) the time evolution of gas tem- 
perature is calculated by solving the radiative transfer equation in- 
cluding the energy deposition of radioactive elements, (ii) the radial 
density profile of various metals is properly considered, and (iii) 



unmixed and uniformly mixed cases in the He core are considered. 
In the unmixed case, the original onion-like structure of elements 
is preserved, and in the mixed case, all the elements are uniformly 
mixed in the helium core. 



It should be mentioned here again that TF01 assumed an adi- 
abatic cooling in the ejecta and adjusted the adiabatic index 7 to 
1.25, referring to the formation episode of dust grains observed in 
SN 1987A. As pointed out bv lKozasa et al] ll989l) . the condensa- 
tion time as well as the resulting average size of dust grains strongly 
depend on the value of 7. SN 1987 A is somehow a peculiar SN in 
the sense that the progenitor is not a red supergiant but blue super- 
giant: see lArnett et alJ Jl989^ for the details. Using SN 1987A as 
a template may not be appropriate in comparisons to supernovae 
as a whole. Therefore in this paper, we use the result of N03 as a 
standard model. 



We should note that N03 also assume the complete forma- 
tion of CO and SiO molecules, neglecting the destruction of those 
molecules, i.e., no carbon-bearing grain condenses in the region of 
C/O < 1 and no Si-bearing grain, except for oxide grains, con- 
denses in the region of Si/O < 1, The formation of CO and SiO 
may be incomplete because of the destruction by energetic electron 
impact within SNe. TF01 treat both formation and destruction of 
CO and SiO, finding that both are mostly destroyed. The decrease 
of CO leads to the formation of carbon grains, which could finally 
be oxidised with available oxygen. The destruction of SiO could 
decrease the formation of grains composed of Si02, MgSi03, and 
Mg2Si04, and increase other oxidised grains and Si grains. Obser- 
vationally, it is still a matter of debate if CO and SiO are efficiently 
destroyed or not. For the detailed discussions on this issue, see Ap- 
pendix B of N03. 
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Figure 2. Evolution of the mass for each grain species. The thick lines indicate the total amount of dust mass. We assumed the closed box model and a 
constant star formation rate. The mass fraction of each grain species is constant in time and is given by the size distribution function (Fig. IT]. 



2.1.2 Dust grain species produced by N03 model 

In the unmixed ejecta, a variety of grain species (Si, Fe, Mg2SiC>4, 
MgSi0 3 , MgO, A1 2 3 , Si0 2 , FeS, and C) condense, and in 
the mixed ejecta, in contrast, only oxide grains (S1O2, MgSi0 3 , 
Mg2Si04, AI2O3, and Fe 3 C>4) form. This is because carbon atoms 
are consumed to form CO molecules. We summarize the species 
formed in SNe in Table Q where the species marked with a circle 
are relevant for unmixed and mixed SNe. The size of the grains 
spans a range of three orders of magnitude, depending on the grain 
species. The size spectrum summed up over all the grain species 
has a very broad distribution, and very roughly speaking, it might 
be approximated by a power law. This size distribution is different 
from that of the SN II calculation of TFO 1 , which has a typical sizes 
of 300 A for amorphous carbon and 10-20 A for oxide grains. 

In this work, we adopt the representative progenitor mass of 
SNe II as 20 Mq . N03 have shown that the size distribution of each 
grain species is almost independent of the progenitor mass, if the 
SN type is fixed (i.e., SN II or PISN). We examine the unmixed and 
mixed cases. The size distributions of dust grains in the mixed and 
unmixed cases calculated by N03 are shown in Figure □ We used 
these size distributions after binning with a bin width of 0.2 dex 
for our calculations. Throughout this work, we assume a uniform 
and spherical grain. It should be mentioned that a different shape 
of dust grains in SN is suggested (e.g., Dwek 2004). 

2.2 Star formation, chemical evolution, and dust production 

For constructing the chemical evolution model of a young galaxy, 
we adopt the following assumptions: 

(i) We use a closed-box model, i.e., we neglect an infall and 
outflow of gas in the scale of a star-forming region. 

(ii ) For the initia l mass function (IMF), we adopt the Salpeter 
IMF <SaipeteJl955t) 

0(m)<xm" 2 - 35 (1) 



Table 2. Mass fraction of dust species (unmixed). 



Species 


Mass fraction 


C 


0.102 


Si 


0.303 


Fe 


0.050 


FeS 


0.156 


AI2O3 


0.002 


MgSiO a 


0.010 


Mg 2 Si0 4 


0.222 


SiO z 


0.099 


MgO 


0.056 



Table 3. Mass fraction of dust species (mixed). 



Species 


Mass fraction 


AI2O3 


0.001 


MgSiOa 


0.105 


Mg 2 Si0 4 


0.202 


Si0 2 


0.583 


Fe 3 4 


0.108 



with mass range of (mi, m u ) = (0.1 Mq, 100 Mq). 

(iii) We neglect the contribution of SNe la and winds from low- 
mass evolved stars to the formation of dust, because we consider 
the timescale younger than 10 9 yr. 

(iv) The interstellar medium is treated as one zone, and the 
growth of dust grains by accretion is neglected. Within the short 
timescale considered here, it can be assumed safely (see, e.g., 
IWhittell992l pp.223-224). 

(v) We also neglect the de struction of dust grains within the 
young age considered (see e.g. Jjones et alll 996). 
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Figure 3. Extinction efficiency Q(a, A) for the grain species produced by the supernova model of Noza waet alj 120031) . We show Q(a, A) for the radius 
a = 10~ 7 , 10 — 6 , 10 — 5 and 10 — 4 cm by the curves from bottom to the top in each panel. 



© 0000 RAS, MNRAS 000, 000-000 



6 T. T. Takeuchi et al. 



Table 4. The specific heat of each species. 



Species 


Debye temperature 


Ref a 


iV Ml-nm 


C(T) model 1 




[K] 




[10 22 cm~ 3 ] 




c 


863 


1 


11.4 


Multi-Debye 


Si 


616 


2 


5.00 


Debye 


Si0 2 




1 


7.96 


Multi-Debye 


Fe 


500 


3 


8.50 


Debye 


FeS 


384 


4 


6.62 


Debye 


Fe 3 4 


511 


5 


28.4 


Debye 


AI2O3 


1030 


6 


23.5 


Debye 


MgO 


762 


7 


10.7 


Debye 


MgSiOs 




1 


9.57 


Multi-Debye 


Mg 2 Si0 4 




1 


9.63 


Multi-Debye 



6 



a References for Debye tempera ture s: CD | Draine <fc l20Ql|). (2) 
[ Qganov. Brodhold. & Pried feOOCl). (3llSiethoff & Ahlbornl <1996l). (4) 
brgnvold et alJ il99lh, (51 IShepherd et alJ <199ll) . f61 iMezzasalmd 
feOOOh . (7l lHama & Suitol<1999l) . 

^ 'Debye' means that the specific heat C(T) is modeled by the clas- 
sical Debye model, while 'Multi- Debye' mean s that the multidimen- 
sional Debye model proposed by Draine & Li 1 2001 ) is adopted. For 
the species to which we adopt the multidimensional Debye model, we 
do not give the classical Debye temperature. 

(vi) We assumed a constant SFR for simplicity. 

using these assumptions, we calculate the chemical evolution. De- 
tails of the formulation are presented in Appendix lAl 

The evolution of the total dust amount for SFR = 1 M© yr -1 
is shown in Figure|2| The dust mass fraction for unmixed and mixed 
cases are summarized in Tables [2] and [5] respectively. Comparing 
the evolution of dust mass given by T03, dust mass starts to accu- 
mulate later than the case of T03, and gradually approaches the T03 
result toward the age of about 10 8 Gyr. This difference is caused by 
the different formul ae we ad opted in the calculations of stellar life- 
time: the formula o f lSchaered l 20021 gives a longer lifetime for the 
same stellar mass than that of Inoue. Hiras hita. & Kamav^ 1:2000) 
which is used in T03. At 10 8 Gyr, T03 and present work yield the 
same dust mass. 



2.3 SED construction 

In this subsection, we present the construction of the SED from 
dust. Since most of the formulations are in parallel to those in T03 
and T04, here we list the essence of our calculation. 

(i) Stochastic heating of very small grains 

Very small grains cannot establish thermal equilibrium with the 
ambient radiation field, which is called stochastic heating (e.g., 
lKriigelll2003l) . To treat this effect, we applied the Debye model 
to the specific heat of the grain species as discussed in T03 and 
T04. We adopt a multidimensional Debye model (e.g., Drai ne & Lil 
l200ll) (DL011 for carbon (C) and silicate (Si0 2 , MgSi0 3 , and 
Mg2SiC>4) grains. For other species, we adopt the classical three- 
dimensional Debye model with a single Debye temperature. The 
specific heat model is summarized in Table[4] 

(ii) Emission 

The emission from dust is calcul ated in the same w ay as 
T03/T04, basically according to lDraine & AndersorJ fl9851. Total 
dust emission is obtained as a superposition of the emission from 



each grain species. The total mass of each grain component is given 
by N03. With this value and material density of each species (Ta- 
ble □ , we can determine the normalization of the dust size distri- 
bution (see AppendixlBl. 

We constructed Q(a,X) of each grain species from available 
experimental data via Mie theory (e.g.. lKriigelll2003l Chap. 2). 
The extinction effici encies are presented in Figure [5] T03 adopted 
Drain e & LeeNl9 841 for the optical properties of these species, be- 
ing different from the present work. On the contrary, we treat de- 
tailed species in the present work. Hence, the mass fraction dis- 
tributed to carbon and silicate are significantly reduced comparing 
to the case of T03. In the IR, Q(a, A) of these grains is much larger 
than other grains at the whole range of grain size. As a result, the 
sum of the contributions of all grains are less than that of the case 
where only two species are considered as T03. 

(iii) Extinction 

Self-absorption in the MIR for a very optically thick case is 
treated by a thin shell approximation, in the same manner as T04. 



3 RESULTS 

3.1 Evolution of infrared SED 

We first show the evolution of the IR SED of forming galaxies 
based on our baseline model. For these calculation we adopted the 
star formation rate SFR = 1 Mq yr" 1 . We adopt rsF = 30 pc 
and 100 pc, the same as those used in T03. These values are rele- 
vant when describing 'dwarf-like' young galaxies. The former suc- 
cessfully reproduce various properties of active star-forming dwarf 
galaxies like SBS 0335—052, whereas the latter is representative of 
the si ze of quiescent dwarf galaxies like I Zw 1 8 1 Hira shita & Hunt! 
l2004ll . We will revisit these representative dwarf starburst galaxies 
in Section|4| 

The results are presented in Figures El and Fi gurel4*lis the 

SED of a galaxy with tsf = 30 pc, while Figure [5] is the one 
for rsF = 100 pc. We calculated the evolution of the SED in 
the age range of 10 6 ' 5 -10 8 yr. In a very young phase (age = 
10 6 75 -10 7 25 yr), unmixed-case SED has an enhanced N-MIR 
continuum. After 10 7 ' 25 yr, the N-MIR continuum is extinguished 
by the self-absorption in the case of rsF = 30 pc. In contrast, the 
self-absorption is not significant for rsF = 100 pc. 

In both cases, the SEDs have their peaks at a wavelength 
A ~ 20-30 /im, which is much shorter than those of dusty giant 
galaxies at 2 = 1-3 detected by SCUBA. At submillimetre wave- 
lengths, the shape of the continuum is very similar to each other 
for mixed and unmixed cases, though the unmixed case predicts 
stronger fluxes by a factor of two. This is explained by the grain 
size distributions of the unmixed and mixed cases: carefully ex- 
amining FigureQ we see that the largest-size dust grains are more 
abundant for the unmixed-case size distribution. It means that larger 
grains contribute more significantly than smaller ones. Since larger 
grains radiate their energy at longer wavelengths, the total SED has 
a stronger FIR-submm continuum. 

As a whole, our model SEDs are less luminous than those of 
T03 by a factor of 2-3. This is due to the difference in consid- 
ered grain species, and more importantly, also due to the improve- 
ment of the treatment of Q(a, A) at the UV; Q(a, A) was approx- 
imated to be unity in T03. Also T03 treated the size dependence 
of Q(a, A) by a simple scaling law, which was not very accurate. 
As seen in Figure|3| the size dependence of Q(a, A) cannot be de- 
scribed by a simple scaling law, especially at the MIR (see, e.g., Si, 
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Figure 4. The evolution of the infrared (IR) spectral energy distribution (SED) of a very young gal axy. The size of the star forming region, rgF, is 30 pc. 
Solid and dotted lines represent the SEDs based on the unmixed and mixed model of dust production i Nozawa et al. 2003), respectively. For Figuresl4landl5l 
a constant star formation rate of SFR = 1 AfQyr - 1 is adopted. 
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Figure 5. Same as Figure|4] but the size of the star forming region rgp is 100 pc. 
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Figure 6. The contribution of each dust species to the total SED of a young galaxy. In this figure we set the galaxy burst age to be t = 10 6 ' 75 yr and the radius 
of the star forming region to be rgp = 30 pc. Left panel is the SED calculated from the unmixed dust production, and right panel is obtained from the mixed 
dust production. 



Fe, and FeS). We, in this work, used the exact value of Q(a, A) in 
all the wavelengths. Since some elements have a very small value 
for Q(a,X) at the UV regime, the total absorption probability of 
UV photons is smaller than that of T03. A strong N-MIR contin- 
uum at a very young phase and the dependence on rsp are qual- 
itatively consistent with the previous result of T03. However, the 
dust grains expected from N03 model predict a peak of the SED 
at shorter wavelengths than those of T03. In addition, a relatively 
small amount of small dust grains also makes the N-MIR SED 
weaker than the T03 one. 



of Fe304 is large, the resultant SED for the mixed case does not 
have a strong N-MIR continuum. 

At A = 10-20 /im, several species contribute equally to 
the SED in unmixed case, while SIO2 play a dominant role for 
mixed case. At longer wavelengths, the main contributors to the 
SED are Mg2SiC>4 and amorphous carbon for the unmixed case, 
and MgSiC>3 and Mg2Si04 in mixed case. In addition, especially 
at A > 200 fim, FeS (unmixed) and FeaCU (mixed) are also im- 
portant because of their shallow slope of the extinction efficiency 
at longer wavelength regime. Note that, for the mixed case, Fe3C>4 
contribute to both NIR and FIR. 



3.2 Contribution of each species to the SED 

The dust species treated here are much more detailed than those 
of T03. In order to see the detailed contributors to the total SEDs, 
we show the individual SED of each species in Figure|6| We show 
the SED with the age of burst 10 6 ' 75 yr, where the effect of self- 
absorption is negligible. 

Both in unmixed and mixed cases, the smallest dust grain 
species is AI2O3 (Figure^. However, because of its small fraction 
in mass (Figure[2j, AI2O3 contributes to the total SED very little. 
The main contributor to the (unextinguished) N-MIR continuum is 
Si for unmixed case and Fe304 in mixed case, respectively. These 
species are the second smallest grains. For the unmixed case, since 
the emissivity of Si grains with the size of ~ 10~ 7 cm is much 
smaller than that of Fe304 at A > 1 ^m, the grain temperature 
becomes much higher than that of Fe3C>4, which makes the very 
strong N-MIR continuum emission. On the contrary, since Q(a, A) 



3.3 Opacity and its evolution 

The extinction curve A\ is obtained as 

A x = 1.086r dust (A) . (2) 

In this work we concentrate on the extinction curve in the IR. De- 
tailed discussions of the extinction curves in the optical wavelength 
regime are given in Hirashita et al. ( 2005|). 2 

FigureQpresents the evolution of the IR opacity for a young 
galaxy with the age of 10 6 ' 75 -10 8 yr (from the bottom to the top). 
For the case of rsp = 30 pc, the extinction becomes almost unity 



2 Since we assume a very simple geometry, the extinction curve is directly 
related to the 'attenuation curve', in which the geometrical effec t is als o 
included (for a thorough discussion on this subject, see, e.g., Inoue 2005). 
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Figure 7. The evolution of the opacity of a very young galaxy at IR wavelengths. Left panel shows the IR extinction curves calculated with the star forming 
region size rgp = 30 pc, and right panel shows the curves with rgp = 100 pc. From the bottom to the top, the ages of galaxies are 

10 6.75 yr> 10 7.0 yr> 

10 7 ' 25 yr, 10 7 ' 5 yr, 10 7 ' 75 yr, and 10 8 yr, respectively. The solid and dashed curves represent the IR extinction curves based on the unmixed and mixed 
model, respectively. 



at 10 fim when the age is 10 7 ' 25 yr. Hence, after this age, the N- 
MIR regime of the SED becomes optically thick, and this makes the 
continuum extinguished in Figure|4| On the other hand, for tsf = 
100 pc, it remains optically thin even at the age of 10 s yr at 10 /im. 

In Figure^ we find an interesting difference between the ex- 
tinction curves of the unmixed and mixed cases. Though their be- 
havior is similar to each other at wavelengths A > 10 /im, there is 
no dip at A < 10 /J,m in the curve of the unmixed case. This makes 
the extinction stronger for the unmixed-case galaxies at N-MIR. In 
addition, the MIR bumps at 10 and 20 /im are less prominent for 
unmixed case. 

Comparing these IR extinction curves with T03 model, we 
find that the amount of the extinction is smaller than that of T03 by 
a factor of 2-10, depending on the wavelength and unmixed/mixed 
production. This is explained by the discussion presented in Sub- 
section [23] T03 supposed silicates and carbon grains as the con- 
stituents of dust grains from SNe. Since their Q(a, A) is larger than 
the other species predicted by N03, the total extinction becomes 
smaller in this work. 

As for the shape of the curve, the extinction curve for the 
mixed model resembles that of T03 extinction curves, while the 
curve for the unmixed model is qualitatively different at MIR, be- 
cause of the weak MIR bump features and the lack of the MIR dip 
mentioned above. 

Figure [3] demonstrates the contribution of each dust grain 
species to the total IR extinction curves. The adopted age is 10 7 yr 
as a representative, but the contribution is constant with time. 

For the extinction curve of the unmixed case, bump features 
at MIR are relatively weak (Figure 0. Similar to those of the 
Galactic IR extinction curve, they are due to silicate grains. In the 



unmixed-model extinction, the bumps are dominated by Mg2SiC>4, 
and the contribution from MgSiOs is not important. We also find a 
strong contribution from amorphous carbon grains which have very 
smooth dependence of Q(a, A) on wavelength. They makes the sili- 
cate features rather weak, because of their comparable contribution 
with Mg2SiO<4. We note that the amorphous carbon does not have 
any dip at A ~ 5-10 /im. Further, large grains of Si have a large 
Q(a, A) at these wavelength regime. We find that both the amor- 
phous carbon and Si grains equally contribute to the total extinction 
at N-MIR, and these species plug up the MIR dip of Mg2SiC>4. At 
NIR wavelengths close to 1 /im, FeS also contributes to the total 
extinction. 

We next see the mixed case. In this case, MgSiC>3 and 
Mg2Si04 equally contribute to the total extinction around the MIR 
bumps. In contrast to the unmixed case, there is no species which 
fill up the dip of silicates, hence we see clear silicate bump features 
and dip in the IR extinction curve. In the N-MIR, main contributors 
to the opacity are Fe304 and Si02- The contribution from AI2O3 
is not important at any wavelength. 



4 DISCUSSION 

4.1 Nearby forming dwarf galaxies 

It is still a difficult task to observe galaxies in their very first phase 
of the SF, especially to detec t their dust emission d irectly. Along 
the line of studies made by Hirashit a et alj 120021) . T03 consid- 
ered two local star-forming dwarf galaxies, SBS 0335—052 and 
I Zw 18, for investigating the dust emission from very young small 
galaxies. Here we revisit these representative star-forming dwarfs 
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Figure 8. The contribution of each dust species to the total IR extinction curve of a young galaxy. In this figure we set the galaxy burst age of t = 10 7 yr. Left 
panel is the extinction curve calculated from the unmixed dust production, and right panel is obtained from the mixed dust production. 



with our framework. Since a recent observation of SBS 0335—052 
by Spitzer has been reported iHouck et alJl2004 . it is timely to 
reconsider these 'textbook objects' with the new data. In addi- 
tion, understanding the SEDs of these objects will shed light to 
the physics of interstellar mat ter and radiation of high-,z ga laxies 
also via empirical stu dies (e.g. jTakeuchi. Yoshikawa. & Ishiil2003t 
iTakeuchi et alJ2005h . 3 



4.1.1 SBS 0335-052 

SBS 0335-0 52 is a local ga laxy (~ 54 Mpc) with SFR = 
1.7 M Q yr -1 jHunt et all200ll) and extremely low metallicity Z = 
1/41 Zq. This galaxy is known to have an unusual IR SED and 
strong flux at N-MIR. It has a very young starburst (age < 5 Myr) 
witho ut significant underlying old stellar population I Vanzi et al. 
2000). T03 have modeled the SED of SBS 0335-052 and reported 
a good agreement with the available observations at that time. 

However, iHouck et all |2004 ) presented new data of the MIR 
SED by Spitzer, and reported a deviation of the model by a factor 
of two or three. Their observation indicated that SBS 0335—052 
has even more FIR-deficient SED than ever thought. Hence, it is 
interesting to examine whether our present model can reproduce 
the extreme SED of this galaxy. We are also interested in the possi- 
bility to determine which of the two pictures, unmixed or mixed, is 
plausible, by a direct measurement of the SED. 

We show the model SEDs for SBS 0335-052 in Figure^ We 

3 We should keep in mind that the following discussion is limited to 
a particular class of dwarf galaxies dominated by newly formed stars. 
For the modeling of the S ED of normal star-forming dwarfs, see, e.g., 
iGalliano et aH<2003tl2005l) . 



have calculated the SED for tsf = 10, 20, and 30 pc both for un- 
mixed and mixed cases. The SFR is fixed to be 1.7 Mq yr -1 , and 
the age is 10 6 ' 5 yr. Solid squares are the observed data from ISO, 
while open squares represent the upper limits obtained from IRAS 
and ISO observations. Open triangles are the expected c ontribution 
of dust emission calculated by the recipe of Ijov & Lesterl ll988l) . 
Filled triangles depict the measured SEP of SBS 0335-052 by 
Spitzer IRS taken from lHouck et alj 12004 in the MIR. The filled 
triangle at FIR is the estimate fro m the radio observ ation of 
iHunt et alJ 120041) . which is used by IHouck et all l2004 . Details 
of the other observational data are found in Section 4 of T03. 

In the FIR regime, ou r model SEDs are c onsistent with the 
strong constraint given by IHouck et alJ l2004h . both for the un- 
mixed and mixed cases. This is because our present model pre- 
dicts a peak of the SED at shorter wavelengths than that of T03. 
At MIR, though we cannot give an excellent fit to the observed 
data, the model SEDs roughly agree with them, and the unmixed- 
case SEDs with tsf = 10-20 pc give a better fit. The very strong 
N-MIR continuum of SBS 0335-052 is well reproduced by the 
SED of the unmixed dust production picture. On the contrary, the 
SEDs of the mixed case seriously underpredicts the observationally 
suggested N-MIR continuum. Since we can also calculate the ex- 
tinction value A\, we can distinguish between a mixed-case SED 
and heavily extinguished unmixed SED without ambiguity. 

In summary, the SED for the unmixed dust production with 
rsF = 10-20 pc yields a reasonable fit to the latest observations 
by Spitzer. This result suggests that we may determine the dust 
production (unmixed or mixed) of SNe through the observation of 
the N-MIR SEDs of forming galaxies. As for SBS 0335-052, the 
unmixed dust production is suggested to take place. In a previ- 
ous work, we have m odeled a high- 2 (z = 6.2) extinction curve 
(Maiolino et al. 2004) and found that it is well represented by dust 
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Figure 9. The model for the SED of a nearby star-forming dwarf galaxy SBS 0335—052. We adopt the star formation rate SFR = 1.7 Mq yr - 1 iHunt et a?] 
l200ll) . The galaxy burst age is found to be t = 10 6,5 yr and the radius of the star forming region is set to be rgp = 10, 20, and 30 pc. Solid squares are the 
observed data from ISO, while open squares represent the upper limits obtained from IRAS and ISO observations. Open triangles are the expected contribution 
of du st emission calculated by the recipe of Joy & Lester 1 1988). Solid triangles are taken from recent observation of SBS 0335—052 by Spitzer iHo uck et alj 
|2004|) . Left panel is the SED calculated from the unmixed dust production, and right panel is obtained from the mixed dust production. Clearly, the mixed dust 
production poorly reproduces the observed SED of SBS 0335—052. 



produced in unmixed SNe jHirashita et ailF2 005 ) . This conclusion 
is strengthened by our result above. 

The dust mass calculated by our model at this age of 
SBS 0335-052 is 1-2 x 10 3 Mp, c onsist ent with the observation- 
ally estimated value by D ale et all 1200 ll) . We note that the mass 
estimation is strongly dependent on the assumed dust species and 
their emissivities, and grain size distribution. As we discussed in 
Section|3| the continuum radiation in the N-MIR is dominated by 
stochastically heated dust emission, which is completely different 
from modified blackbody. Therefore, when we try to estimate the 
dust mass, we must take care to determine the corresponding grain 
properties, i.e., radiative processes and grain which are related to 
the observed SED of galaxies. 

4.1.2 IZwl8 

I Zw 18 is the most metal-deficient star-forming galaxy in the lo- 
cal universe (Z = 1/50 Zq). We adopt its distance of 12.6 Mpc 
16stlirl200rj) . and the SFR of 0.04 M Q yr" 1 . The existence of the 
underlying old stellar population is suggested for this galaxy, but 
their contribution is not dominant lHunt etall2003h . Thus, it is not 
unreasonable to adopt our present model to predict the IR SED of 

IZw 18. 

ICannon et alj 1200 2) estimated the total dust mass Md us t = 
2-3 x 10 3 M e from HST WFPC2 narrow-band imaging. In 
our model, the dust mass reaches this value at burst age of 
10 7 ' 25 -10 7 ' 5 yr. This age is consistent with th e obser vationally sug- 
gested age of its major SF (e.g.. IHunt et alJl2003h . ICannon et alJ 
(2002) obtained the extinction of Av = 0.5 mag for some patches 
in this galaxy, which is reproduced by our model with the above 
age. 

We show our model prediction for the SED of I Zw 1 8 in Fig- 
urefTol We present the SED for the age 10 7 ' 25 yr and 10 7 ' 5 yr. The 



NIR flux measurements are taken from lHunt etafl 120031) . I Zw 18 
is not detected at all the four IRAS bands, hence we calculated the 
upper limits in the M-FIR from the sensitivity limits of IRAS ob- 
servation. At 850 /xm, the u pper limit is converted from the data 
reported bv lHunt et all 120051) . 

The IR luminosity is slightly smaller than the previous predic- 
tion of T03. We also find that the FIR peak wavelength is located at 
A ~ 60 /im, shorter than that of the T03 model SED. In addition, 
the peak intensity is larger for the mixed model, since the mixed 
model produce about 30 % larger amount of dust than the unmixed 
model does (see Fig. [2j. Since the system still remains optically 
thin at this age because of the large tsf, this difference is directly 
reflected to the difference of the peak intensity. 

The low opacity of I Zw 18 also leads to a low flux at N-MIR 
in the model SED in Figure fTol both for unmixed and mixed cases 
compared with the T03 model SED. As mentioned by T03, the NIR 
continuum of I Zw 18 is expected to be dominated by stellar and 
nebular continuum radiation, and the contribution of dust may be 
negligible. 

We expect that these features will be found commonly in 
young galaxies if they are optically thin in the IR. Hence, it may 
be useful to obtain some suggestions for the observational strat- 
egy of young galaxies at high-z. We will discuss a well known real 
sample of such population of galaxies, Lyman-break galaxies, in 
the next subsection. 



4.2 Lyman-break galaxies 

Lyman-break galaxies (LBGs) are one of the m ost well-s tudied cat- 
egories of high-z star-forming galaxies (e.g., Steid eTet alJI 1999t 
2003). Even in LBGs, there is clear evidence that they contain 
non-negligible amount of dust (e.g., lAdelberger & Steidelll200ol 
ICalzettilboOll) . A high dust temperature (> 70 K) is suggested 
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Figure 10. A prediction for the SED of a nearby 'quiescent' star-forming dwarf galaxy I Zw 18. We adopt SFR = 0.04 M0yr _1 and rgp = 100 pc. The 
age of the star formation is set to be 10 7,25 yr (dashed lines) and 10 7 ' 5 yr (solid lines). Since the NIR continuum of I Zw 18 is expected to be dominated 
by stellar and nebular continuum radiation, the contribution of dust may be negligible, and it is not necessary to reproduce these NIR observed fluxes by the 
present model. 



by subsequen t studies lOuchi et alJ 1 19991 IChapman et ail [2000 ; 
ICalzettfcOO il : Is awickil200 ll) . 

In order to make a consistent picture of the dust emission from 
LBGs, T04 applied the T03 model and made various predictions for 
these galaxies. They also considered a power-law dust size distri- 
bution to examine its effect on the resulting SEDs. The T04 model 
reproduced the known observational properties of the dust emission 
from LBGs. 

In this paper, we investigate the expected appearance of the 
LBGs with an improved dust grain formation of N03. 

4.2.1 Evolution of the SED of LBGs 

We set the input parameters of the SED model for LBGs as 
follows. The SFR of LBGs spreads over the range of SFR ~ 
1-300 Mpyr" 1 with a median of SFR ~ 20M Q yr _1 (e.g., 
lErb et alJl2003h . A constant SFR up to the age of 10 9 yr is sug- 
gested to be a good approximation (e.g., Baker et al. 2004). Thus, 
the basic framework of the T03 model is also valid for LBGs. In 
this work, we consider the moderate case of SFR = 3OM0yr _1 
over the age of 10 6 ' 5 -10 8 yr. 

The most important information to calculate the IR SED is 
the effective size of the star forming region, but it is the most un- 
certain quantity (see discussion of T04). Since the mean half-light 
radius of LBGs is estimated to be ~ 1.6 kpc from HST observa- 
tions <Erbetall2 003). we use the galaxy radius as the radius of a 
star-forming region, and set rsF = 2 kpc according to T04. 

The SED evolution of LBGs is shown in Figure fTTI Note that 
the scale in the ordinate of Figure ITU is different from those in 
Figures UJ and |5] Aside from luminosity, the behavior of the SED 
evolution of LBG is qualitatively very similar to that of a dwarf- 
like forming galaxy with rsF = 100 pc (i.e., a weak NIR con- 
tinuum and a strong M-FIR emission is commonly seen in both). 
This is because both the UV photon number density and dust num- 
ber density are in the same order for an LBG and dwarf galaxy 



(rsF = 100 pc). Since the IR luminosity depends on the dust col- 
umn density, the total luminosity is different, but the shape of the 
SED is determined by the balance between the densities of UV pho- 
tons and dust grains, they are similar with each other as for the 
shape. 

The difference between the SEDs of the unmixed and mixed 
cases is small, and they give almost the same result especially in 
the M-FIR. This is explained as follows: the most prominent differ- 
ence between the two is the N-MIR continuum, which is produced 
mainly by stochastic heating of grains. If the UV photon density 
is low, they are almost always in the lowest temperature state be- 
cause the incidence of a UV photon is less frequent. In this case the 
N-MIR continuum from stochastically heated dust does not con- 
tribute the continuum significantly, and the global SED shape is 
determined by larger grains in the equilibrium with the ambient 
UV radiation field. Under such condition, they are similar in shape 
with each other. Thus, we see that it will be difficult to determine 
which scenario of dust production is correct from the observation 
of LBGs. 



4.2.2 Observability of the dust emission from LBG 

The flux density of a source at observed frequency, u h B is obtained 
by 

= (1 + z)L (1+z>obB = (1 + z)L„ cm 

° ba A7Yd L (z) 2 47Td L (2) 2 ' K ' 

where d^(z) is the luminosity distance corresponding to a redshift 

2, and ^obs and v arn are observed and emitted frequency, respec- 
tively. 

We show the observed IR/submm SEDs of LBGs at z = 2, 

3, and 4 in Figure fT2l For simplicity we only show the SED with 
the age of 10 s yr. The thick black short horizontal lines indicate the 
3-cr detection limits for 8-hour observation by ALMA (Atacama 
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Figure 11. The evolution of the SED of a typical Lyman-break galaxy (LBG). The SFR is set to be SFR = 30 MQyr - x , and the size of the star-forming 
region rgp = 2 kpc. 




Figure 12. The expected flux densities of an LBG. In this figure we set the galaxy burst age of t = 10 8 yr. Confusion limits of Herschel and detection limits 
of ALMA 8-hour survey are also shown by gray and black thick horizontal lines, respectively. 



Large Millimeter Array). 4 Here we assumed 64 antennas and three 
wavelength bands, 450, 850, and 1080 fim. We also show the 3-cr 
source confusion limit of Herschel 5 at 75, 160, 250, and 350 /im 
bands by thick gray horizontal lines. These limits are based on 'the 



4 URL: |http : / /www . alma . inf o/| 

5 URL: |http : / / www .rssd.esa. int /herschel /| 



photometric criterion' of jLg^a che et all2003l) . See also llshii et alJ 
i2002h : lTakeuchi & Ishii <2004al) . 

As discussed in T04, the detectability of LBGs is not strongly 
dependent on their redshifts. The predicted IR luminosity is factor 
of two smaller than that of T04, hence it is more difficult to de- 
tect LBGs than the discussion in T04: we need roughly four times 
longer integration time than that derived from T04, if we fix all 
the other conditions. Detection at 350 fim seems impossible for 
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moderate-SFR LBGs. However at longer wavelengths, if the age 
> 10 8 yr and SFR > 30 M yr~ 1 , LBGs can be detected at a wide 
range of redshifts in the submm by ALMA deep survey. In the FIR, 
Herschel will detect the dust emission from LBGs at z ~ 2, but dif- 
ficult at higher-z. Since the SED does not differ for the unmixed or 
mixed dust formation, our evaluation of the observability of LBGs 
holds both for the mixed and unmixed dust formation. 



4.3 Toward higher redshifts 

Based on the above discussions, we give a brief consideration on 
the observation of very high-z galaxies (z > 5) here. Direct ob- 
servation of such galaxies are of vital importance to explore the 
physics of galaxy formation. 

First, assume a LBG at these redshifts. Figure [T3l shows such 
a situation. We assumed the same physical parameters for LBGs. 
From Figure fT3l we find that LBG-like objects can be detectable 
even at z ~ 10 if SFR > 30 M Q yr _1 . At z ~ 20, it may be more 
difficult, but might be possible to detect LBGs with very high SFR, 
by an ALMA ultra-deep survey with much longer integration time 
than eight hour. 

However, it may not be reasonable to assume a galaxy like 
LBGs at z ~ 10 because they are rather massive system in general 
(M star — 10 10 Mq), and, from the modern cosmological view- 
point, such a massive system is very rare at such a high-z, young 
universe. 

Here we discuss an object less massive than LBGs. In modern 
hierarchical structure formation scenarios, it would be more reason- 
able to assume a small, subgalactic clump as a first forming galaxy. 
Consider a dark halo of mass ~ 10 9 Mq, then it is expected to con- 
tain a gas with mass ~ 10 8 Mq. For this purpose, we calculate the 
SEDs for a dwarf star-forming galaxy. If gas collapses on the free- 
fall timescale with an efficiency of esf (we ass ume esF = 0-1), 
we ob tain the following evaluation of the SFR Hira shita & HunJ 
12004 : 

n H 
100 cm" 3 



SFR-O.lfg^f M ^ )( WH -V /2 [Meyr- 1 ] , (4) 

V o.i y ^io 7 Mq) Vioocm- 3 ; 1 oy J ' 



and 

where nu is the hydrogen number density. Then we have 
SFR * 0.1 (^) ( JfeY /2 (y^) " 72 [Moyr- 1 ] . (6) 



Vo.iy ^io 7 M e 



If we consider M gas ~ 10 s Mq, we have SFR ~ 
3(r SF /100 pc)" 3/2 MQyr" 1 . In addition, as we see below, an 
extremely hi gh-z galaxy observ ed by HST has a very compact 
morphology iKneib et alJ 12004 1. We also mention that, from a 
theoretical side, high-z galaxies are sug gested to be dense and 
compact compared to ne arby galaxies jNorman & Spaan3ll997t 
Hiras hita & Ferraral2002l) . Hence, it may be reasonable to assume 
a same type of galaxy as local dwarfs discussed in Section 14.2.1 1 
Thus, we consider a dwarf galaxy with SFR = 10 Afgyr - as an 
example, and we adopt rsF = 30 pc and 100 pc. The age is set to 
be 10 7 yr. If the age is older, they will become easier to detect if 
a constant SFR takes place. We show the expected SEDs for such 
galaxies at z = 5, 10, and 20 in Figures l"i"4l and IT5l As expected, 
it seems almost impossible to detect such objects by Herschel or 
ALMA. 



There is, however, a hope to observe such a small forming 
galactic clump directly: gravitational lensing works very well as a 
natural huge telescope. For example, MS 1512— cB58 is an idea l 
case of a lensed LBG l Nakani slfietai]|l997t iBaker et alj|2004) . 
and more recently, a small star-forming ga laxy has been detected 
by H ST and Spitzer observations iKneib et alJl2004j: lEgami et alJ 
2005). Hence, we can expect a forming galaxy which is intrinsi- 
cally small and faint but magnified by a gravitational potential of a 
cluster of galaxies. If we assume a lens magnification factor of 40, 
such a small galaxy becomes detectable. This is depicted by the thin 
lines in Figures fT4l and fT5l Since the expected SED of such a com- 
pact dwarf galaxy has a strong MIR continuum at their rest frame, 
it can be feasible to detect at the FIR in the observed frame. Indeed, 
at z = 5 they appear above the confusion limits of Herschel. Since 
the confusion limit of Herschel may be, in fact, very difficult to 
reach in actual observations, a cooled FIR space telescope is more 
suitable for such observation, and this will be a strong scientific 
motivation for a future project like SPICA. 6 In the FIR wavelength 
regime (50-200 /im) , the SEDs of unmixed and mixed cases might 
be distinguishable. If a FIR spectroscopy is performed, which pic- 
ture of dust production is plausible might be examined. At higher- 
z, they can be detectable by ALMA survey. Even at z ~ 20, they 
can be detected by a standard 8-hour survey of ALMA, if a lens- 
ing takes place. It is interesting to note that, at z ~ 10, the peak 
of their SED locates at the wavelength range where the detection 
limits of both Herschel and ALMA are not very good. A new ef- 
fective (probably space) facility should be developed to overcome 
this difficulty in observation. 

For this line of study, we must estimate how frequently such 
lensing events occur for high-z objects. Suppose a cluster of galax- 
ies at z\ ~ 0.1-0.2 whose dynamical mass M dyn is 5 x 10 14 M Q 
and whose mass distribution obeys the singular isothermal sphere. 
We denote the strong lensing cross section, i.e., the area of the re- 
gion in the source plane for which the resulting magnification by 
a cluster is larger than /i, as <j(> u). IPerrotta et alJ 120021) pre- 
sented er(> fj) as a function of Md yn for zi cns = 1.0. Since cr(> 
fi) oc Di s 2 [Dis is the angular-diam et er distance between the lens 
and t he source (e.g.. lAsadalll998l) l ICovone. Sereno. & de Ritisl 
2005), we can convert their result to our condition and obtain 
u(> 10) ~ 30 arcsec 2 on the source plane. 7 This result is almost 
independent of the source redshifts. Setting the limiting flux density 
S v = 1 fj,]y and using the number counts of Hirashi ta & Ferrara 
( 2002) for galaxies at z > 5, we have an expected number of galax- 
ies suffering a strong lensing to be ~ 1-3. Thus, we expect at least 
a few strongly lensed IR galaxies to this survey depth. 

To have a concrete idea for the observational strategy, we 
must consider an important aspect of the lensing. It is known 
that the gravitational lensing affects the galaxy number counts in 
two opposite directions: i) it magnifies the flux density of galax- 
ies, and ii) it stretches the sky area and decreases their num- 
ber density on the sky. Thi s is ca ll ed the 'ma gnification bias' 
iFort Mellier. & Dantel-Fortl Il997t iBezecourt. Pello. & Soucaill 
1998). The unlensed number count Nq(> S v ) is generally approx- 
imated by a power law 



N (> S v ) oc Sj 



(7) 



Then, by denoting the magnification factor as fi, we obtain the 



6 URL: |http : / /www, ir . is as ■ jaxa ■ jp/ SPICA/ index . html| 

7 On the image plane, the area corresponds to ~ Aic(> /i). 
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5.0 SFR= 30 Moyr-' 
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Figure 13. The expected flux densities of an LBGs hypothetically located at z = 5, 10, and 20. The assumed parameters are the same as those for Figure lb?! 
Solid lines represent the SEDs for the unmixed case of dust production, and dotted lines for the mixed case. We put the same observational detection limits as 
Figure fT2lfor Herschel and ALMA. 
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Figure 14. The expected flux densities of a dwarf star-forming galaxy located at z = 5, 10, and 20. Solid lines represent the SEDs for the unmixed case of 
dust production, and dotted lines for the mixed case. We assumed SFR = 10 M0yr — 1 , and rgp = 30 pc. Thick lines are the expected SEDs for an unlensed 
galaxy, while the thin lines are the ones for a gravitationally lensed galaxy by a magnification factor of 40. Again the detection limits are the same as those of 
Figure fT2l 
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Figure 15. The same as Figure [l4l except that rgp = 100 pc. 
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lensed number counts iVi cns (> S v ) as 

AW(> &) = N {> &)/u" (/3+1) (8) 

(e.g., iFort. Mellier. & Dantel-Fortl Il997t 

Takeuchi. Yoshi kawa^^^^neharal200ol) ! Equation shows 
that whether the number counts increases or decreases depends 
on the slope of the counts, /3: if /3 < — 1 lensing works as 
an enhancement of the counts. Since the slope of the number 
count s of such galaxies i s expected to be steeper than —1 (e.g., 
Hiras hita & Ferrard fooCOl the source density on the image plane 
will increase (see also Perrotta et al. 2002). Thus, we are confident 
that the lensing will work as a really useful tool to detect forming 
small galaxies with dust emission. 



parameters of LBGs, the unmixed and mixed picture does not af- 
fect the appearance of the SED. Hence, if we combine our model 
with the knowledge obtained from optical observations, the observ- 
ability of LBGs at submillimetre wavelengths are robust conclusion 
independent of the details of the SNe dust production theory. 

Finally, we considered the observations of forming galaxies 
at 2 > 5. If there exist LBG-like galaxies at these redshifts, they 
can be detected at z < 10 by ALMA 8-hour survey if SFR > 
30 M©yr~ . For small forming galaxies or subgalactic clumps, it 
is almost impossible to detect their intrinsic flux by ALMA or by 
Herschel. However, the gravitational lensing is found to be a very 
effective tool to detect such small star-forming galaxies at z > 5. 
If we consider a compact dwarf galaxy with SFR ~ 10 A/oyr -1 
and rsF — 30 pc, we can expect a few strongly magnified galaxies 
behind a typical cluster of galaxies at z ~ 0.1-0.2. 



5 CONCLUSION 

Dust plays various important rol es even in the v ery early phase of 
galaxy evolution (e.g.. lHirashita & Ferraral2002h . In such a young 
phase with a typical age less than 10 9 yr, dust is predominantly 
supplied by SNe. With the aid of a new physical model of dust 
production by SNe developed by N03, we constructed a model of 
dust emission from a very young galaxies according to T03. 

N03 carefully took into account the radial density profile and 
the temperature evolution in the calculation of the dust formation in 
the ejecta of SNe II and PISNe. They also showed that the produced 
dust species depends strongly on the mixing within SNe. We treated 
both unmixed and mixed cases and calculated the IR SED of young 
galaxies for both cases. 

The SEDs constructed from N03 dust production are less lu- 
minous than those by T03 model by a factor of 2-3. This difference 
is due to the improvement in the treatment of Q(a, A) at UV and the 
considered grain species. The SED for the unmixed case is found 
to have a strong N-MIR continuum radiation in its early phase of 
the evolution (age < 10 7 ' 25 yr) compared with that for the mixed 
case. The N-MIR continuum is due to the emission from Si grains, 
which only exist in the species of the unmixed dust production. 

We also calculated the IR extinction curves for young galax- 
ies. N03 dust gives a weaker extinction than that of T03 model be- 
cause of the small relative number of very small dust grains. This 
is also the cause of the smaller IR luminosity of the present model. 
For the unmixed case, NIR extinction is dominated by large grain of 
Si and amorphous carbon, and silicate features are less prominent 
compared to the curve given by T03. On the contrary, the extinction 
curve of the mixed case has a similar shape with that of T03. 

The SED of a local starbursting dwarf galaxy, SBS 0335—052, 
was calculated. Recent Spitzer observations iHouck etal]|2004 
have implied a hotter dust than previously thought. Our present 
model SED naturally reproduces the strong N-MIR continuum and 
the lack of cold FIR emission of SBS 0335-052. We found that 
only the SED of unmixed case can reproduce the N-MIR contin- 
uum of this galaxy. Hence, as for SBS 0335—052, the unmixed dust 
production is preferred. It will be interesting to proceed this line of 
study for higher- z galaxies. 

A prediction for the SED of another typical nearby star- 
forming dwarf galaxy, I Zw 18, was then made. Since I Zw 18 is 
dominated by very young star formation, we may adopt the model 
to this galaxy. Using the dust species from N03 with this model, we 
find a weaker FIR emission than that of T03. The N-MIR contin- 
uum is also expected to be much weaker than that of T03 SED. 

We also calculated the evolution of the SED of LBGs. For the 



ACKNOWLEDGMENTS 

First we are grateful to the anonymous referee for many insightful 
comments which have improved the clarity of this paper much. We 
deeply thank Akio K. Inoue, Jean-Paul Kneib and Giovanni Covone 
for illuminating discussions. We made extensive use of the NASA 
Astrophysics Data System. TTT, TTI, and HH have been supported 
by the Japan Society of the Promotion of Science. TK is supported 
by a Grant-in-Aid for Scientific Research from JSPS (16340051). 



APPENDIX A: DETAILS OF THE CHEMICAL 
EVOLUTION 

Here we present a detailed description of the chemical evolution 
we used in this work. Full treatment of the present chemical evolu- 
tion model including the feedback, dust destruction, and more com- 
plex SF history is given elsewhere (Nozawa et al., in preparation). 
Within the framework of the closed-box model, the total baryonic 
mass of a galaxy, Mt, is conserved as 

M T = M gas (0) = M gBS (t) + Mdust(t) + AW(t) + AW(t) , (Al) 

where M gas (t) is the gas mass, Md us t (t) is the dust mass, M s tar (i) 
is the stellar mass, and M rem (t) is the mass of stellar remnants in 
the galaxy at an age of t. We start the calculation from homoge- 
neous pristine gas, i.e., M dust (0) = Mstar (0) = M rom (0) = 0. 
The time evolution of the mass of interstellar medium, Mism(£) = 
M gas (t) + Mdust (t) is given by 

— — — = —ip(t) + / [m - m rcm (m)] ip{t - r m )<j>(m)dm , (A2) 

* J mi 

where tp(t) is the SFR at an age t, <j>( m ) is the IMF normalized 
to unity in the mass interval [mi, m u ] (subscripts 1 and u mean 
lower and upper mass limit, respectively), m le m(m) is the mass of 
a single stellar remnant resulting from a progenitor star with mass 
m, and r m is the lifetime of a star with mass m. 

For the IMF, as we mentioned in the main text, we adopt the 
Salpeter IMF with mass range of (m h m u ) = (0.1 M Q , 100 M ). 
The lifetime r m is evaluated by a fittin g formula o f zero-metallicity 
tracks with no mass loss, given by Schaerer ( 2002). For SNe II (pro- 
genitor mass 8-40 Mq), we approximate the remnant mass m rcm 
by the l inear fit based on the num erical result of SN explosions 
given by Umeda & Nomoto ( 2002): 

m Iem (m) = 0.06m + 0.93 (8 M < m < 40 M Q ) . (A3) 
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For other progenitor mass ranges, we set m rem (m) = at m < 
8 Mq (no remnants), and m lem (m) — m at 40 Mq < m < 
100 Mq (all the mass is swallowed by a black hole). In this work, 
we do not treat PISNe in the calculation. By Equation <A3> . the 
time evolution of M tem (t) is expressed by 



dM r , 



dt 



x (m)i/](t — T m )(j}(m)dm . 



The time evolution of A/dust (t) is, then, obtained as 



dM d 



dt 



M duBt (t) i 
'MisM(i) 1 



tp(t) + / m Aust (m)ip(t - 



(A4) 



z )(f)(m)dm 



where md us t(m) denotes the dust mass produced by a progeni- 
tor with mass m. We apply the dust formation model of N03 for 
Tidust- As we mentioned in § i|2,ll we used the result of progenitor 
mass 20 Mq for the fraction and the size distribution of dust grain 
species, hence the mass fraction and size distribution are constant 
in time along with the evolution. 



APPENDIX B: NORMALIZATION OF THE NUMBER 
DENSITY OF GRAIN SPECIES 

In this Appendix, we show how to obtain the normalization of the 
dust size distribution dNi/dai, 



^^-da = Afi(a)da 



(Bl) 



where subscript i denotes the species of dust, and fi is the mass 
fraction of dust of species i, given by N03, and A is the normaliza- 
tion constant. We denote the material density of i-species dust by 
pi, and the SF region radius as rsF- The total mass of dust species 
i, Mi, is then written as 



M, = Pi 



Ana d dNi 
3 da 



da — pi 



Afi(a)da 



(B2) 



The total dust mass of all the species, M, is 

[it A y" 



M = y> 4ttA^ f /i(a)a 3 dfl = ^y^ p . I /i(a ) a 3 da . ( B3) 



Hence 
A = - 



3M 



47r > pi I fi(a)a 3 da 



(B4) 



Thus we can determine A from the available data. Note that only 
M depends on time t, hence the time dependence of A is given 
only through M. The total number of grains, Ni, is expressed as 

V, = / ^da t . (B5) 
dai 
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